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We introduce an approach to quantum information processing where the information is stored 
in the motional degrees of freedom of nanomechanical devices. The qubits of our approach are 
formed by the two lowest energy levels of mechanical resonators which are tuned to be strongly 
anharmonic by suitable electrostatic fields. Single qubit rotations are conducted by radio frequency 
voltage pulses that are applied to individual resonators. Two qubit entangling gates in turn are 
implemented via a coupling of two qubits to a common optical resonance of a high finesse cavity. 
We find that gate fidelities exceeding 99% can be achieved for realistic experimental parameters. 
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Mechanical oscillators are among the most elementary 
structures that are studied in physics. Nonetheless they 
have properties that make them very useful for techno- 
logical applications. Their vibrational modes can for ex- 
ample undergo millions of oscillations before the motion 
is eventually damped and they can couple to electromag- 
netic fields in a very broad frequency range via their po- 
larizability. Whereas the latter property has prompted 
significant effort to build optical to microwave frequency 
converters [TJ [5], the enormous Q-factors of mechani- 
cal oscillators have been exploited to demonstrate ap- 
proaches to mechanical computers in the classical domain 
[2l |4] . The fascinating properties of nanomechanical os- 
cillators furthermore motivated intense research activity 
towards exploring their quantum regime which has very 
recently lead to breakthroughs in cooling such oscillators 
to their quantum ground states jSH?]. 

Here we introduce an approach to quantum informa- 
tion processing with mechanical degrees of freedom by 
making use of the aforementioned exquisite properties of 
nanomechanical oscillators. The device we envision, con- 
sist of an array of doubly clamped nanobeams that all 
couple to a common resonance mode of a high finesse op- 
tical cavity, see figure [T] for an illustration and possible 
setup. Each nanobeam can furthermore be manipulated 
individually with electrostatic and radio frequency fields 
that are applied via small tip electrodes. A suitable setup 
for an implementation are for example carbon nanotubes 
that couple to the evanescent field of a whispering gallery 
mode cavity [SHU], c.f. figure [l] 

In our scheme, the constant electric fields applied to 
individual nanobeams make the mechanical spectrum of 
each beam anharmonic and allow to tune the respective 
transition frequencies. For large enough fields, the an- 
harmonicity of the mechanical spectrum becomes compa- 
rable to the line-width of the optical cavity resonances. 
By driving the cavity with a coherent input that is ap- 
propriately detuned in frequency from the closest cavity 
resonance, one can ensure that only one transition be- 
tween eigenstates of the mechanical motion couples to 
the cavity. Here, we choose the laser drives such that the 
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FIG. 1: Top row: Illustration of nanomechanical qubits in- 
teracting via a common photon mode: a) The deflection of 
one resonator locally changes the energy density of the pho- 
ton field and thus causes a force onto other resonators, b) By 
properly tuning the qubit frequencies, noninteracting subsets 
can be defined. Bottom row: Possible implementation of the 
device we envision, c) The mechanical vibration of doubly 
clamped carbon nanotubes couples optomechanically to the 
evanescent field of a whispering gallery mode in a high finesse 
micro-toroid. d) Electrostatic and radio frequency fields can 
be applied to each nanotube (black) individually via tip elec- 
trodes (blue). 



cavity couples to the transition between the ground and 
first excited state of the nanobeam. These two states, 
denoted |0) and |1), form our nanomechanical qubit. 

As compared to previous approaches |15l 116] . our 
scheme uses the intrinsic nonlinearity of nanomechanical 
resonators which allows to selectively address the qubit 
transitions and results in very high gate fidelities. Lo- 
cal operations respectively single qubit gates are imple- 
mented by applying suitable voltage pulses via the tip 
electrodes. The optical cavity mode in turn mediates 
interactions between any pair of mechanical qubits. If 
such an interaction is active for a suitable time range 
and combined with pertinent local operations it can im- 
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plement a fundamental entangling gate, e.g. the so called 
iSWAP gate [19]. For many qubits, one can selectively 
apply such an iSWAP gate to any desired pair of qubits 
by tuning them to a transition frequency ujq while the 
remaining qubits are tuned to a markedly different transi- 
tion frequency cos via suitable voltages at the respective 
tip electrodes. As we show in detail below, at the end 
of the iSWAP gate on the selected qubits all remaining 
qubits can be made to return to their initial state and 
effectively undergo an identity operation. After success- 
fully completing a quantum algorithm, the state of each 
qubit can be read out by tuning individual qubits to dis- 
tinct transition frequencies and performing spectroscopy 
with a weak probe laser, c.f. [13 ] 114 ) . 

With all the above ingredients, the device we envi- 
sion satisfies the requirements for implementing quantum 
computing 20 : 1) It has well defined qubits formed by 
the two lowest energy eigenlevels of strongly anharmonic 
nanomechanical oscillators and is scalable since multiple 
nanobeams can couple to one high finesse optical mode. 
2) It is initializable in the state |0, 0, . . . , 0) by cooling the 
mechanical motion to the ground states, e.g. via side- 
band cooling [6, 7 . 3) The extremely high Q-factors of 
nanomechanical oscillators together with a cryogenic en- 
vironment ensure that the coherence times of the qubits 
greatly exceed gate operation times. 4) A universal set 
of gates can be implemented. 5) Spectroscopy on indi- 
vidual qubits together with single qubit rotations allows 
for high quantum efficiency, qubit-specific readout. 

The system We consider a system consisting of N 
nanobeams that are clamped at both ends and thus 
feature an intrinsic nonlinearity that originates in the 
stretching of the material associated with its deflection 
[2"Tl |2"2"] . In terms of phonon creation (annihilation) op- 
erators bj (bj), the Hamiltonian of one nanobeam reads, 

ff£°i = &^&}& i + ^(&} + & J )\ (i) 

where w m is the resonance frequency of the harmonic mo- 
tion and A = is^zPM the nonlinear contribution. Here 
ft depends exclusively on the dimensions and material 
properties of the beam whereas xzpm denotes the ampli- 
tude of its zero point motion. All nanobeams are subject 
to static and radio frequency electric potentials gener- 
ated by tip electrodes in the close vicinity of each beam 
[T31 IT41 IT71 [P5] and couple to a common optical mode 
of a high finesse cavity. The Hamiltonian of the entire 
electro-opto-mechanical setup in a frame that rotates at 
the frequency of the drive laser reads, 

H = - HAa^a + hg^X ^ X 3 + hgX c ^ X j 
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where A is the detuning between driving laser and cav- 
ity resonance and Xj = (bj + bj)/y/2 the deflection of 
beam j. The photon operators a and a? have been 
shifted by their steady state values, a —¥ a + a and a 
negligible term cx a'a^j^j has been dropped. Hence 
X c — (a* a + aa' )/V2\a\ is a photon quadrature with 
a = 2A+i K i where Q is the drive amplitude of the laser 
and k the photon decay rate of the cavity. 

The potentials V* v ' z (t) = V^ z + V*J' z (t) describe 
constant, V^' z , and time dependent, Vjf' z (t), gradient 
forces caused by the voltages applied to the tip elec- 
trodes. The constant parts tune the equilibrium posi- 
tions of the mechanical oscillators via V^q and their spec- 
trum via VjQ j whereas the time dependent parts Vj 
can implement single qubit rotations. We always choose 
VjQ = ~h~9~^ such that the nanobeams remain unde- 
flected. g = 2|a|a;zpMGo is the optomechanical coupling 
that can be controlled by the amplitude of the laser drive, 
where Go = duj/dX is the cavity's frequency shift per 
resonator deflection. We use the radio frequency fields 
Vj t \(t) and the coupling g as controls to perform gate 
operations. 

In a realistic experimental situation, the mechanical 
motion will be subject to damping at a rate 7, n and cavity 
photons will be lost at a rate k. The full dynamics of our 
system that takes these incoherent processes into account 
can thus be described by the master equation, 

Q = -i[H,Q] (3) 

3 

where n is the thermal phonon number at the environ- 
ment temperature T and the dissipators read 2?j,(j/) = 
2?/£>2/ t ~ y f yg - gy f y and V^(y) = 2y^gy - yy^ g - qyy^ 

Nanomechanical qubits The qubits we consider are 
formed by the two lowest energy levels of our nanome- 
chanical beams. Mechanical resonators feature a small 
instrinsic anharmonicity in any deflection mode that can 
be enhanced by electrostatic gradient forces, see [13l [14] 
for details. Suitably tuned electrostatic fields generate an 
inverted harmonic potential of the form V cx —X 2 that 
counteracts the harmonic part of the elastic restoring 
force. The electrostatic potential thus reduces the stiff- 
ness respectively softens the mechanical resonance mode 
and hence reduces its frequency. Since the deflection per 

— 1/2 

phonon, xzpm, is proportional to w m a reduction of 
w m causes a significant enhancement of the nonlinearity 
as A oc x| PM oc lu~ 2 , c.f. [13] [14] . 

It is useful to consider a "tuned" mechanical Hamil- 
tonian, H m j = H^j + J^j VfflXj, that includes the 
constant parts of the V? . We write H m j and the de- 
flection Xj of each mechanical resonator in the eigen- 
basis of H m j, i.e. H m ,j = ^2 n E n \ri)j{n\j and Xj = 
J2nm X rim\n)j(m\j, where H m<j \n) = E n \n). E n and 
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FIG. 2: a) Nonlinear phonon spectrum in the presence of a 
softening field Vq . For Vq > Y^mfii f ne nanobeam enters 
the bistable buckling regime, b) Effective nonlinearity of the 
phonon spectrum S21—S10. A sufficiently high value of S21— <5io 
suppresses unwanted processes of the type |11...) — > |20...}. 

X nm can be found numerically, see FIG. [2^i. Here we 
assumed all resonators to be identical and thus skipped 
the subscript j in the E n and X nm . A generalization 
to inhomogeneous arrays is straight forward. Provided 
the nonlinearity per phonon is large enough and the me- 
chanical motion is cooled to the groundstate, we can re- 
strict our analysis to the two lowest energy levels |0) 
and |1) which form our qubit. The nonlinearity of the 
qubits, i.e. the fact that the qubit transition energy 
5 10 = Ei — Eq differs significantly from transition en- 
ergies between higher states S nm = E n — E m , especially 
from S21 (see FIG. [2]b) , ensures that the restriction to 
these two levels remains a very accurate approximation 
throughout the gate sequences. To explain the gate op- 
erations in detail we now switch to an interaction picture 
with respect to Hq = — HAaJa + . H m j. 

Local Operations and Single Qubit Gates Local oper- 
ations, that is rotations about the o~ x -, a y - and ovaxes 
are realized by applying time dependent gradient forces 
as encoded in the potentials V^ v (t) and V?(t) in eq. 

cr z -rotations are obtained by temporarily shifting the 
qubit transition frequency <5io which makes the qubit ro- 
tate at a different rate and hence collect a phase shift. 
This is achieved by tuning the softening fields that con- 
trol the mechanical frequency to a different value for a 
suitable time which adds a term V Z X 2 to the Hamilto- 
nian that shifts the qubit frequency #10 — > <5io + 8$ (t) 
(We skip the qubit index j throughout the discussion of 
local operations). For J S^dt — </>/2, this procedure im- 
plements the operation 

e -i/v-(t)x'd*/n „ e -i^/2 =_[^. (4) 

a x - and <jj,-rotations can be implemented by apply- 
ing gradient forces related to V xy . If the potential 
V xy (t) is a pulse that is modulated by an oscillation 
at the qubit frequency V xy (t) = cos(S w t + 9)V xy ', 

one finds V xy cos(5 w t + 6)e iHot / h Xe~ iHa ^ h V xy x 
x^l [co S (0) ( a oi +a Wj +isin (0) (3.01 _ o-io-j] 5 where 

a ab = \a)(b\ and we have dropped off-resonant contri- 
butions. Thus, for 9 = and J dtV xy (t) = (f>/X 01 one 



finds 

e -i/dtv-»(t}x RJe -^/a s _ [ ^ ]Tj (5) 
and for 9 — ir/2 and j dtV xy (t) = 4>/Xq\ 

e -i/dtV-*(t)JC fae -ur v */2=_ [ ^ ]iri (g) 

Effective interactions and two- qubit gates In order to 
implement entangling two-qubit operations we employ an 
indirect interaction between mechanical resonators that 
is mediated by an optical cavity mode [53]. Here, the 
laser is sufficiently far detuned from any resonance, such 
that g <C I A ± 5 nm \. For an initial state in the qubit 
subspace formed by |0) and |1), this condition ensures 
that the dynamics is restricted to this subspace. To per- 
form a gate operation on two qubits, e.g. with indices 
j = {1,2} ("gate qubits"), we tune those two qubits to 
a common transition frequency ujq, while tuning the re- 
maining qubits ("saved qubits") to a sufficiently differ- 
ent transition frequency cjg such that \ujq — cjg| 3> g 
and interactions between gate qubits and saved qubits 
are strongly suppressed. To explain the working princi- 
ple of this gate, we consider a scenario where V^f' z = 0. 
An adiabatic eliminiation of the photons together with a 
rotating wave approximation yields the effective Hamil- 
tonian 

H c s ~ Hq (t) + Hs ■ (7) 

with H G = -H0^ (a? 1 af + H.c.) - Y%=i H G,i and 
Hs = -fc^E^aW 1 ^ +H.c.) -E i>2 Hs,i, 
where H G/8>i = ft£ (c° /s af + c^a} 1 ) and of = 

\a){b\i. Here, ATq/s denotes the displacement matrix ele- 

x 2 

ment X01 and c G ^ s = Y^, m A +7" * ne interaction induced 
phase shift for the gate qubits and saved qubits, respec- 
tively. This phase shift can be reversed after the gate 
operation by appending a suitable ovrotation. 

The time evolution of the Hamiltonian Q can be 
used to perform an iSWAP operation on distinct gate 
qubits, while the other qubits are unaffected. This can 
be achieved by decomposing the iSWAP operation into 
V ±iSWAP operations, using the identify 
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Here, the V±iSWAP operations are implemented by 
choosing laser drive pulses such that A ^°a / g(t) 2 dt = 

±^ respectively. During the V iS WAP operation on the 
gate qubits, the saved qubits are subject to the dynamics 
generated by H savc (t) (compare FIG.[T]d). Yet since there 
are no local — [tt}z operations on the saved qubits, this 
evolution is reversed during the \/ — iSWAP operation on 
the gate qubits and the saved qubits return back to their 
initial states. 



Numerical results We analyzed the fidelites of the 
quantum gates we propose by numerically solving ^ in- 
volving two and four qubits where we included the low- 
est three levels |0), |1), |2) for each resonator respectively 
qubit. To find an estimate for the fidelity of the gate op- 
erations we compute the fidelity f(a, p) = Tr(-\/ \fp V\fp) 
of the desired target state a with the actual state p that 
results from the full dynamical evolution for several ini- 
tial states. In order to confirm that relative phases of 
the states evolve as desired, we use all states of the 
form \</>) = (\j,k,...) + |0,0,...))A/2 (j, k, . . . = 0,1) 
as initial states and average the fidelity / over these, 
F = f(< J ,p)\ aU \ c f i y FIG. jij shows the resulting gate error 
£ = 1 — F as a function of the gate time Tq, the laser 
detuning A and the mechanical and optical dissipation 
rates 7 m and k. To demonstrate the scalability of our 
approach we present results for two qubits (solid lines) 
as well as for four qubits (red dots). For the four qubit 
case, the gate is applied to two qubits while the remain- 
ing two return to their initial state. These results cleary 
show the excellent scaling properties of our approach. 

The result given in FIG. H are found for (10, 0)- 
carbon nanotubes of length L = 300 nm and radius 
R = 0.39 nm, coupled to the evanescent field of a mi- 
cro toroid cavity, see [13j [14] for details. By soft- 
ening the mechanical resonance frequency of the gate 
qubits down to uj g /2tt = 26.6 MHz, we achieve a non- 
linearity of (S 2 i - S w )/2n = 2.71MHz. The optome- 
chanical coupling of carbon nanotubes can be dramati- 
cally enhanced by employing a cavity resonance in the 
vicinity of recently demonstrated excitonic resonances 
[30t [31] . We choose a cavity frequency that is far enough 
from such a resonance to sufficiently suppress additional 
absorption and estimate that an optomechanical cou- 
pling rate of g/2n = 8.73 MHz can be achieved with a 
laser drive of 5.30 W input power that is detuned by 
A/2-7T = 0.399 GHz from the cavity frequency. For a 
gate time of Tq = 1.49 us such coupling suffices for the 
ViSWAP operation asuming a rectangular pulse. Fur- 
thermore we asume an ambient temperature of 20 mK, a 
mechanical Q-factor of 5 • 10 6 for the tuned resonator [TO] 
and a total photon decay rate (including intrinsic cav- 
ity losses and losses induced by the exciton resonance) 
of k/2tt = 133 kHz [21 [25]. Note that for our con- 
ditions with k 2 /A 2 = 8 x 10 -8 this corresponds to an 
absorbed power of ~ 18 nW that is compfortly compat- 
ible with a cryogenic environment |26j . Moreover even 
higher finesses and lower absorption could be attained 
in crystalline resonators [27]. For tuning the mechan- 
ical resonances, electrostatic fields of 76 V/um are re- 
quired, whereas for the local operations (T x ,y\<p\, we find 
radio frequency field amplitudes of ~ 1 V/pm, depend- 
ing on operation time and angle (f>. For the static fields 
that compensate the photon induced shift of the equilib- 
rium position during an entangling operation, 6.1 V/pm 
are needed and the cr z [<j)] operations require a gate time 
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FIG. 3: Error £ — 1 — F of gate operations. Except for 
the quantities on the horizontal axis all parameters are as 
in the setting described in the main text. Solid lines show 
results for two and red dots for four qubits. Highlighted red 
dots correspond to the parameter example discussed in the 
text, a) Error as a function of the gate time for local and an 
entangling operations. Note that changing the gate time also 
changes the required coupling rate, b) Error as a function 
of the laser detuning, c) and d) illustrate the influence of 
damping for the photons and the resonators, repectively. 



> 4>/u>q to avoid buckling the resonator. 

Conclusions In summary, we have introduced an ap- 
proach to quantum information processing with nanome- 
chanical qubits. Our approach is realizable with current 
or near future experimental technology. Importantly the 
nanomcchanical resonators in our approach show very 
promising scalability properties and are not vulnerable to 
fluctuations of background charges such as superconduct- 
ing qubits. Alternative platforms for an implementation 
could consist of stiff nanobeams with high optomechani- 
cal couplings such as photonic crystal nanobeam cavities 
in diamond |28j . For building even larger scale devices 
our scheme could also be integrated into optomechanical 
networks of multiple cavities [29] . 
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SUPPLEMENTARY INFORMATION 

The optomechanical coupling of the nanotube to the 
cavity is mediated by gradient forces, originating from 
the evanescent cavity field and acting on the nanores- 
onator. Therefore, the strength of optomechanical cou- 
pling is proportional to the polarizability a of the nan- 
otube. Here, we make use of the fact, that at the scale 
of |iin optical wavelength, the polarizability can be sig- 
nificantly enhanced due to the presence of excitonic reso- 
nances [3D] . Near an excitonic resonance of frequency tu c 
and line width T c , the polarizability can be estimated by 



a(ui) 



c 2 f 



1 



iwIY 



(9) 



where / « 10 — 100 is a typical oscillator strength in these 
systems [32]. With fiw 8 al,4 eV and 1/T e w 1 ns [3DJ, 
this leads to a potential amplification of the polarizability 
and hence, optomechanical coupling strength, by a fac- 
tor > 10 6 (at the maximum of 5J(a)) as compared to the 
static value of 10 — 100 A . However, to minimize ad- 
ditional photon losses introduced by the imaginary part 
3(a), it is preferable to choose the cavity to be detuned 
from the exciton resonance (compare Fig. [4|. We choose 
a detuning lu c — cu c = 10~ 4 w o , which yields an enhance- 
ment factor of roughly 2 x 10 4 , and additional losses of 
k c R3 60 kHz. These losses add to the intrinsic cavity 
losses to yield a broadened cavity linewidth k, which we 
considered in our calculatios. 
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FIG. 4: Sketch of the real and imaginary part of the dy- 
namical polarizability at an excitonic resonance. We choose 
a cavity frequency lu c that is sufficiently detuned from the 
resonance to suppress losses induced by Q(a). 



